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We present an extended study of anyonic Luttinger liquids wires jointing at a single point. The 
model on the full line is solved with bosonization and the junction of an arbitrary number of wires is 
treated imposing boundary conditions that preserve exact solvability in the bosonic language. This 
allows us to reach, in the low momentum regime, some of the critical fixed points found with the 
electronic boundary conditions. The stability of all the fixed points is discussed. 



I. INTRODUCTION 

In the last few years there has been a boom in the 
study of transport properties at the junction of multiple 
quantum wires 1-26 . This interest is largely motivated by 
the fact that junctions of three or more wires would nat- 
urally appear in any quantum circuit. Different frame- 
works have been developed to tackle this complicated 
problem that shows a rich phase diagram. In fact, de- 
spite of the universality in the bulk of the wires, that are 
described by a Luttinger liquid 27 , different conditions at 
the junctions can lead to exotic phase diagrams (as e.g. 
those in Refs. 1,9,23) whose degree of universality is 
not yet understood. According to the Renormalization 
Group (RG) theory of critical phenomena, the low en- 
ergy properties of a gapless system are captured by the 
stable fixed point of the RG flow, independently of micro- 
scopic (non-universal) details of the real system. In view 
of the universality it is worthy to investigate very simple 
models, even exactly solvable, that can have (because 
of symmetry reasons) the same fixed points of the real 
systems. For bulk one-dimensional (ID) models and in 
the case of a single boundary, conformal field theory pro- 
vides a complete classification of the universality classes 
(see e.g. Ref. 28), whose analogous for junctions (or 
star-graph) is not yet known. For all these reasons, we 
investigate in this paper the Tomonaga-Luttinger (TL) 
model on a junction with an arbitrary number n of arms 
as depicted in Fig. 1 (a junction with two wires n = 2 
can be seen as a defect on the line, a problem that has 
been largely investigated 29-40 in the past). To solve this 
problem, at the junction we impose conditions that are 
probably not obvious for an electronic problem, but they 
show the advantage to be exactly solvable. The natural 
hope is that the electronic model, at least for some values 
of the couplings, would be in the domain of attraction of 
the fixed points found here. 

Furthermore we calculate the transport for particles 
with generalized anyonic statistics 41 . The reason for this 
generalization is twofold. On one hand the study of 
ID anyonic model is attracting a renewed interest 42-59 , 
mainly motivated by possible experiments with cold 
atoms 60 . On the other hand, the transport of wires 
joined with a quantum Hall island is driven by anyonic 
excitations 12 . Also in this case we can wonder whether 



the different problems have some common fixed points. 
In ID, anyonic statistics are described in terms of fields 
that at different points [x\ ^ x 2 ) satisfy the commuta- 
tion relations 

&(t,xi)tf(t,x 2 ) = e mK < Xl2) ^{t,x 2 )^{t, Xl ), (1) 

where e(x) is the sign function [e(z) = — e(— z) = 1 for 
z > and e(0) = 0] and x\ 2 — x\ — x 2 . n is called 
statistical parameter and equals for bosons and 1 for 
fermions. Other values of k give rise to general anyonic 
statistics "interpolating" between the two familiar ones. 

The TL model emerges naturally in the description 
of spinless fermions in ID (and so electrons when the 
spin degrees of freedom are not important, but spin is 
also easily introduced in the formalism). In fact, starting 
from fermions hopping on a chain, linearizing the dis- 
persion relation close to the Fermi surface at ±kp and 
taking the continuum limit, one arrives to the standard 
TL Hamiltonian 27 

H = J dx [v F (ij)lid x il)i - Tp* 2 id x ip 2 ) + g+p\ + g-pl] , 

(2) 

where ipi^ 2 {t,x) are the two complex fields representing 
free- fermions left and right movers, Vf is the Fermi ve- 
locity, i.e. the speed of the non interacting fermions, and 

P± (t, x) = [ipl (t, x)tpi (t, x) ± Va (t, x)ip 2 (t, x)] , (3) 

are the two independent charge densities. All the inter- 
action is encoded in the coupling constants g± [often the 
couplings <?2,4 = 2(g+ T <?-) are used]. Eventual irrele- 
vant coupling terms of degree greater than four have been 
dropped. For g+ > g_ the model is repulsive and it is 
attractive in the opposite case. 

A similar reasoning can be repeated for anyonic de- 
grees of freedom and the Hamiltonian is always given by 
Eq. (2), but with ip a satisfying the commutation rela- 
tions (1), ipi with k and tp 2 with —k. Thus, when k = 1 
the model is the well-known fermionic TL model, while 
the bosonic limit k — ► is not well defined in this for- 
malism as will be clearer in the following. We stress that 
this anyonic model is different from the gases discussed 
elsewhere 42,43 ' 46,49 ' 56 , that also have a Luttinger liquid 
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FIG. 1: (Color online) A quantum junction of n wires con- 
nected via a scattering matrix S. 



description. As in the fermionic case, the model is natu- 
rally solved exactly through bosonization 27 . 

This Hamiltonian defines completely the system on 
each wire. To complete the description of the junction 
like the one shown in Fig. 1 we have to define the inter- 
action between the n wires. From an electronic point of 
view it is natural to have a term of the form 1 



(4) 



where a, (3 = 1,2 and i,j = l...n. The matrix B defines 
the boundary interaction among the fields ip. Although 
very natural, this boundary condition is quite compli- 
cated after bosonization, because it involves exponential 
boundary interactions of the bosonic fields. As a conse- 
quence the theory with this interaction term is no longer 
exactly solvable with bosonization, and very smart and 
complicated methods must be employed to extract the 
low-energy behavior from it 1 ' 9 . In this paper we take 
an alternative approach that is to modify the junction 
couplings in such a way to preserve the exact solvabil- 
ity after bosonization 20-22 . The main idea is to impose 
the boundary condition directly on the bosonic degrees 
of freedom trying to have the same symmetries as in the 
■0 counterpart. The two problems can obviously have a 
different structure of fixed points, but, as stressed above, 
the natural hope is that the junction defined by Eq. (4) 
shares some of the anyonic fixed points with the ones 
found here, as it is well known to happen for fermions. 
The clear advantage of our approach is that keeping ex- 
actly solvability, the results are obtained with a relative 
little effort, compared to analogous ones for Eq. (4). 

The paper is organized as follows. In Sec. II we intro- 
duce the anyonic TL model and solve it on the full line. 
In Sec. Ill after introducing the general features of the 



junction and the importance of conservation laws, we first 
present the standard solution on half-line and then gen- 
eralize it to the generic junction. In Sec. IV we study the 
stability of the found fixed points following the RG flow. 
Finally in Sec. V we draw our conclusion and discuss is- 
sues that need further investigation. In two Appendices 
A and B we report the technicality of bosonization and 
the description of the fixed points of the junction. 



II. THE ANYONIC TOMONAGA-LUTTINGER 
MODEL 

As already mentioned, the main goal of this paper is 
to investigate the Tomonaga-Luttinger model on half in- 
finite quantum wires jointing in a single junction. How- 
ever, in order to fix the notations and some basic tools, it 
is instructive to sketch first the solution of model on the 
line. In doing that we will focus on the general anyonic 
solution, which contains the more familiar fermionic one 
as a special case. The model is defined by the Hamilto- 
nian (2) in which the space variable x is integrated on 
the full real axis. The corresponding equations of motion 
are 

i(d t - v F d x )ipi(t,x) = 

2 9+ P+ (t, x)tpt (t,x) + 2<?_ p_ (t, x)tpi (t,x), 
i(d t +v F d x )ip 2 (t,x) = 

2.9+ P+ (t, x)tp2 (*, x) + 2g_ p_{t,x)ip 2 {%x) . (5) 

Bosonization 27 is the basic tool to quantize and solve 
these equations of motion. In fact, the solution can be 
expressed in terms of the right and left-moving scalar 
fields <pr,l- The standard details of the solution can be 
found in textbooks 27 and are reported in appendix A to 
make this paper self contained. The method is based on 
the change of variable 

lpl(t,x) (X : e ^l°VR(vt-x)+T VL (vt+x)] . ^ ^ 

1p 2 (t,x) (X : e ^\TVR.{vt-x)+<"PL{vt+x)] ;j ( ? ) 

where the proportionality constants are explicitly given 
in appendix A, and : • • • : denotes the normal product 
relative to the creation and annihilation operators of tp 
fields, a, t and v are three real parameters to be de- 
termined inserting these expressions in the equations of 
motion. Without loss of generality we take a > and 
assume that 

<7^±T. (8) 

The charge densities take the very simple form 
-1 



p±(t,x) 



[(dip R )(vt-x)±(dip L )(vt + x)]. 

(9) 



2 v ^F(r ± a 
Imposing the current conservation 

d t p±(t, x) - vd x j±(t, x) = , 



(10) 
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one gets the currents 

(dfR)(vt -x)t {d(p L ){vt + x) 



3± (*, x) 



2V^(r ± a) 



(11) 



Using the exchange properties of ip, one can easily show 
that that the field ip\ satisfies the anyonic commutation 
relations given in Eq. (1) with statistical parameter 



2 2 

K = t -a 



(12) 



According to Eq. (8) n ^ 0, that shows explicitly that 
the bosonic limit is not well defined in this context. The 
exchange relations of ip 2 follow from Eq. (1) with the 
substitution n — k, implying that ip a are both anyon 
fields, which become canonical fermions for k = 1. 

The quantum equations of motion are obtained from 
Eq. (5) by replacing p ± (t, x)ip a (t, x) i — p ± ip a :, giving 



t(v — vf)tt z 
a(v + vf)tt 



9-\ 



9- 



t + a t — a 
9+ 9- 



t + a t — a 



(13) 
(14) 



which combined with Eq. (12) determine a, t and the 
velocity v in terms of the coupling constants g± and the 
statistical parameter k. In terms of the variables £± = 
t ± a, one obtains the system of equations 



C+C- = K, 

,2 2 

<+ = v f k + -g+ , 

7T 



Vfk H — 9- 

7T 



with solution 



2 / TTKVp + 2g+ \ 



±1/2 



^(ttkvf + 2g^){-Knvp + 2g + ) 
tt\k\ 



(15) 
(16) 

(17) 

(18) 
(19) 



The relations (18) and (19) are the anyonic realization 
of the well known result valid for canonical fermions in 
the TL model (the traditionally used parameter K in 
our notation coincides for n = 1 with = (+ 2 , for 
comparison in Refs. 9,23 the notation is g — K^ 1 ). The 
stability conditions of the model is 2g± > —itkvf that 
ensures a, r and v to be real and finite. 

From the previously given mapping it is easy to write 
the Hamiltonian in terms of the bosonic fields, obtaining 



h = - I dx [(d x ef + {d xV f} , (20) 



where 



<P{t,x) = - [<p R (vt-x) + ip L (vt + x)] , (21) 
0(t, x) = \ [ip R (vt -x)-<p L (vt + x)] , (22) 



where 9 is the so-called dual field. Notice that the Hamil- 
tonian is slightly different from the usual one in the liter- 
ature because we adsorb the coupling constant g (or K) 
in the definition of the fields. 

It is worth commenting at this point the internal sym- 
metries of the TL Hamiltonian, because they will charac- 
terize the quantization on the junction. The TL Hamil- 
tonian (2) is left invariant by the two independent U(l) 
phase transformations usually denoted as U(l) <g> t/(l): 



i> 



-i(-l)' 



c- 



- e i( - i)a >: 



(23) 
(24) 



In the bosonic language they correspond to the indepen- 
dent shift-invariance of the (compactified) fields <pu.L- 
We will see that on the junction, the left and right movers 
are not independent anymore and the two U(l) symme- 
tries cannot be conserved simultaneously 

One of the main advantages of bosonization is that af- 
ter having solved the equations of motion, it is straight- 
forward to obtain all the correlation functions (also at 
finite temperature) just by commuting the fields ip in the 
exponential forms for ip, using Eq. (A15). In fact, in 
terms of the basic correlator 



V{x) 



1 



i(x — ie) 



(25) 



the zero-temperature (Fock representation) field correla- 
tion functions are 

(il>l{ti,xi)ip 1 (t 2 ,x 2 )) 

= ^-{V(vt 12 -x 12 )} a2 [D(vt 12 +x 12 )Y 2 , 
(ip 2 (h,xi)ip 2 (t 2 ,x 2 )} 

= ^-[V(vt 12 - x 12 )f \v(vt 12 + x 12 )f , (26) 

with Xi 2 — xi — x 2 and t\ 2 =t\ — t 2 . Scale invariance is 
manifest and one can read the dimension of ip a 



dv^ = \{<? 2 +T 2 ) = \(Cl+£). 



(27) 



All the other two-point field correlation functions vanish 
because of Eq. (8) and the neutrality condition (£7(1) ® 
[/(l)-symmetry). Analogously for the [/(l)-density one 
finds 

(p+(h,x 1 )p + (t 2 ,x 2 )) 
= 72^j2 [[ V (vt 12 - x 12 )} 2 + [D(vt 12 + x 12 )] 2 } , (28) 

and straightforwardly the ones for p_ and j± are ob- 
tained. We notice that all these correlation functions 
correctly agree with the general expression for an har- 
monic anyonic fluid 48 with only one harmonic term given 
by the Luttinger mode. 
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The generalization to finite temperature (3 1 (Gibbs 
representation) is simply obtained with the replacement 
V(x) -» T> (x) with 



V {x) 



i(3 I TTX 

sinb I — - K 



(29) 



and introducing the chemical potentials, explicitly 

(lpl(tl,Xl)ll>l(t2,X2))0 

_ _j_ e iM R CT ( t ' t i2-2;i2)+iM L r(wti2+xi 2 ) 

2tt 

x [©^ («t 12 - a; 12 )] CT2 [V p (vt 12 + x 12 )Y 2 , (30) 

and similarly for the other correlations. The right and 
left chemical potentials are 

A*h = ~i i V L =-- + --, (31) 

C.+ s- <■>+ <■>- 

where it and /2 are the ones associated with the (7(1) <8> 
J7(l)-charges. 



III. THE JUNCTION OF 
TOMONAGA-LUTTINGER LIQUIDS 

A. Boundary conditions and symmetries 

After the previous preliminary considerations on the 
line, we investigate below the TL model at a junction 
like the one shown in Fig. 1. In mathematical physics 
literature these junctions are usually called star graphs 
and they represent the building blocks for more general 
"quantum graph" networks (see for a review Ref. 61). 
We now fix all the notation on the junction that we call 
T. We indicate the jointing point of the junction as V. 
Each point P in the bulk T\V (i.e. of the wires) can be 
parametrized by the pair (a;, i), where i = 1, . . . , n labels 
the edge Ei and x € (0, oo) is the distance of P from the 
vertex V along that edge. We stress that, as physically 
suggested, the embedding of T and the relative position 
of the edges in the "ambient space" are irrelevant. 

The dynamics of each wire (edge) is still given by the 
Hamiltonian (2), but now ip a — ip a (t,x,i) and x > 0. As 
already discussed in the introduction, in order to fix the 
solution one must impose some boundary conditions at 
the vertex V at x — 0. The simplest boundary condition 
one can imagine is linear in ip a and is generated by the 
boundary term in Eq. (4) that makes the model non- 
exactly solvable for general couplings (see e.g. Ref. 1 for 
free fermions and also 9 for infinite repulsive coupling). 

An alternative which preserves the exact solvability 
after bosonization has been proposed 20-22 . The main 
idea is to impose the boundary condition directly on the 
bosonic degrees of freedom, selecting those of them which 
ensure unitary time evolution of the fields tp. This is guar- 
anteed only if the boundary conditions are linear in the 



fields tp and its first derivatives. So we can parametrize 
these boundary conditions by a generic n x n unitary 
matrix U 20 < 21 - 62 - 63 



- U)y tp(t, 0, j) - i(I + U)y (d x tp)(t, 0, j)\ 



0. 



(32) 

and A > is a parameter with dimension of mass 
needed to recover the correct physical dimensions. Since 
bosonization expresses physical charges linearly in tp, we 
shall see below that these boundary conditions simply 
state how the charges are parcelled out among the wires 
at the vertex. 

The analysis of the fixed point is greatly simplified if 
we assume time-reversal invariance. This implies that 
the matrix U must be real, that together with unitarity 
leads to a symmetric matrix U, i.e. 



U* =U, 



(33) 



giving a further constraint on the possible boundary 
terms. A non trivial magnetic flux (breaking time- 
reversal) has been considered 9 and resulted in a more 
complicated fixed point structure. When dealing with 
anyon excitation, it would be more natural to consider 
non time-reversal models, because the magnetic field 
needed to produce the anyons breaks the symmetry. 
However this would complicate the analysis and in some 
regime it could be only an irrelevant perturbation. Thus 
in the following we will always assume time-reversal in- 
variance and leave the study of the effect of its breaking 
to a future work. 

The boundary condition (32) is equivalent 20 ' 21 to an 
interaction with a point-like defect localized at the vertex 
of the graph. The scattering matrix associated with this 
interaction i s 20 ' 21 > 62 

S(k) - -[A(I-U) + /^I+U)]- 1 [A(I-U)-/c(I+TLJ)], (34) 

and has transparent physical meaning: the diagonal el- 
ement §u(k) represents the reflection amplitude on the 
edge Ei , whereas (fc) with i =/= j equals the transmis- 
sion amplitude from Ei to Ej . Eq. (34) makes also clear 
the meaning of the boundary terms A and U: for A ^ 
we have §(fc = A) = U, i.e. A fixes the momentum scale 
at which the scattering matrix is given exactly by U. 
By construction the scattering matrix (34) is unitary 

[S(fc)r = [S(fe)]- 1 ) (35) 

and satisfies Hcrmitian analyticity 

[S(ft)]* = S(-ft) . (36) 

Moreover, time reversal invariance (33) implies 

[S(fc)]* = S(fc) . (37) 

For simplicity we assume in this paper that U is such 
that 



/°° Ah 



x > , (38) 
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which guarantees that S(fc) has no bound states (see 64 
for an extension to bound states). 

The boundary conditions strongly influence the sym- 
metry content on the junction. Each symmetry in the 
bulk gives a conserved charge Q [with density p(x, t)] be- 
cause of the Noether theorem. If we want to keep the 
conservation of Q at the junction we must impose from 
the beginning that the currents j(x,t) corresponding to 
the given density p(x, t) are conserved at the vertex. This 
results in i(0> = ® f° r au times. This is the Kirch- 
hoff 's rule, which must be imposed in the vertex in order 
to generate a time-independent charge from a given cur- 
rent. A basic example is given by the energy, that is 
a conserved quantity in the bulk. Because of unitarity, 
the matrix U in Eq. (32) parametrizes all boundary con- 
ditions which ensure the Kirchhoff rule for the energy- 
momentum tensor of ip and thus the time-independence 
of the relative Hamiltonian. This means that there is no 
dissipation at the junction: if the energy flows out from 
one wire should flow in another one. We stress that the 
Kirchhoff 's rule for gapless models on a graph is the gen- 
eralization of the celebrated result that scale invariancc 
implies holomorphic and antiholomorphic components of 
the energy tensor to be equal in boundary conformal field 
theory 28 - 65 . 

Energy is not the only conserved quantity. In our 
formalism it is conserved by construction, but all other 
conservation laws we want to keep on the junction must 
be imposed by hand with appropriate Kirchhoff's rules. 
However it may happen that different conserved currents 
can generate contradictory Kirchhoff's rules, resulting in 
obstructions for lifting all symmetries on the line to sym- 
metries on r 66 . In this case one can preserve on T one of 
the corresponding symmetries, but not all of them. This 
is actually the case for the U(l) ® [/(l)-symmetry of the 
TL model. In fact, the relative Kirchhoff rules generate 21 
the following further constraints on U 



E u(t, o, i) = o ^ e m*o - E u ^ = 1 > ( 39 ) 

i—l i—l i—l 

n n n 

Ei-M,*) = o EM*0 = E u ^ = - 1 - ( 4 °) 



i=i 



i=l 



i=l 



which cannot be satisfied simultaneously. U(l) is linked 
to the electric charge conservation and it is then natural 
to require the conservation of Eq. (39), while breaking 
Eq. (40). However also the opposite prescription has 
some interest. Notice that the duality transformation 
(A12) on T maps the matrix U [and so S(fc, A)] in — U 
[— S(fc -1 , A -1 )]. Consequently duality maps the vertex 
conservation of U(l) in U(l). 

The matrix conductance G of the junction can be 
obtained in linear response theory. Since it involves 
only currents, the calculation is the same as for free 
bosons 21 ' 22 , but with the renormalized current in Eq. 



(11), leading to an overall normalization: 

G=-l T (I-S)=G line (I-S). (41) 
2ttQ + 

Thus the dependence of the conductance on the anyonic 
parameter is only through the renormalization constant 
(+ in Eq. (18). Because of unitarity \Su\ < 1, we have 



< Gjj < 2Gii no . 



(42) 



In the following, we will call conductance G the diagonal 
element Gu in the case it does not depend on the wire 
index i. 

It is worth mentioning that a similar approach (called 
Delayed Evaluation of Boundary Condition) working also 
with fermion boundary conditions has been developed by 
Chamon et al. 9,23 . It basically amounts to leave in the 
half-line, right and left movers unconstrained in the bulk, 
constructing then the tunneling operators, and only later 
choosing an i?-matrix (R for reflection, it can be easily 
rewritten as an S-matrix) such that one of these processes 
pins the correct boundary conditions. In the appendix A 
of Ref. 23 the conductance is written in terms of an n x n 
R, which agrees with the results here and elsewhere 21 ' 22 . 



B. The half-line 

It is instructive to start with the well-known case n = 
1, namely the half-line, since some features of the generic 
junction are already manifest in this case. The matrices 
U and § are just numbers U and S. Setting U = e~ 2lQ , 
we get 



with 



S(k) 



n = Atan(a) , 



k — irj 
k + \rj 



TT TT 

— < a < - 
2 ~ ~ 2 



(43) 



(44) 



As expected the S-matrix (43) corresponds to full reflec- 
tion and describes the mixed (Robin) boundary condition 



{d x <p){t,0)-ri<p(t,0)=0. 



(45) 



The condition (38) implies 77 > or equivalently < 
a < it/2, a — and a = tt/2 correspond to Neumann 
and Dirichlet boundary conditions respectively. These 
two points define the only bosonic scale invariant bound- 
ary conditions on the half-line. Instead of imposing the 
condition (45), we can add a term to the Hamiltonian in 
such a way to generate it as a further equation of motion. 
The resulting total Hamiltonian is 



H Tot = W + r^ 2 (i,0) 



(46) 



with H the bulk term given by Eq. (20), obviously 
defined only on the half-line, i.e. the integral is over 
x G (0, 00). 
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The main effect of the boundary in x = is to couple 
right and left movers by means of the boundary condition 
(45). In particular, at criticality, Eq. (45) implies that 

MO = MO, v = o, (47) 

MO = -MO, V = oo, (48) 

which is the familiar "unfolded picture" 1 for Neumann 
and Dirichlet boundary conditions. The boundary con- 
ditions then forces non-zero mixed commutation relation 
[from Eqs. (A19) and (A20)] between right and left 
movers 

[<Pr(Ci).Vl(6)] = 

-ie(Ci2), V = 0, 

ie(£i 2 ), 7? = oo, (49) 

ie(a 2 )-4i0(ei 2 )e-"« 12 , < 77 < oo, 

while the left-left and right-right ones are the same as in 
the full line. Note that in the right-left commutators it 
appears £12 = vt\ 2 — X12, involving, as expected, the sum 
of distances from the boundary X12 = x\ + X2- 

Although right and left modes are no longer decoupled, 
we can still perform the bosonization program and solve 
the TL model on the half-line. The anyonic exchange 
relations (1) are still valid defining tp a as in Eqs. (6) 
and (7) [but with normalization constants depending on 
the boundary conditions, see Eq. (A18)]. ip a fulfills the 
quantum equations of motion of the TL model restricted 
to the half-line x > 0, with a, t and v given by the same 
expressions (18) and (19) found for the full line. In fact, 
all the local bulk relations of the TL model on the full 
line still hold on the half-line. This will remain true in 
the more general case of a junction made of any number 
of wires. 

The charge and current densities Eqs. (9) and (11) are 
still locally conserved [i.e. Eq. (10) holds for x ^ 0], and 
p± generate the [/(1)®Z7(1) infinitesimal transformations 
(A16). After bosonization, the boundary condition (45) 
can be recasted in terms of physical currents 



j+(*,0) = 0, 
i-(i,0) = 0, 
d x j-(t, 0)-w_(t,0) 



0, 



77 = 0, 

7] = OO, 

< 77 < 00 . 



(50) 



Consequently, the main physical difference between half 
and full line concerns the global charges Q and Q asso- 
ciated to charge densities p + and p_ respectively. The 
boundary spoils the simultaneous conservation of both 
charges, allowing just one linear combination to survive. 
For instance, at the critical point 77 = 0, the boundary 
condition (45) is simply the Kirchhoff 's rule associated to 
the U(l) transformation (23), enforcing the charge den- 
sity current j+ to vanish at the vertex while j_ does not 

j+(i,0) = 0, j_(i,0)^0, for 77-O. (51) 

In this case Q is time- independent, while Q depends on 
time due to a nontrivial charge flow through the bound- 
ary. The critical point 77 = 00 has an opposite behavior, 



preserving the U(l) transformation (24), and breaking 
(23). For generic finite 77 > 0, it is easy to see that 
the U(l) symmetry is always conserved while [7(1 ) is 
broken 34 . As already pointed out, this symmetry break- 
ing from U(l) ® [7(1) to a subgroup U(l) is a general 
unavoidable feature of junctions of any number of wires. 

This boundary symmetry breaking is even more visi- 
ble in the correlation functions. In addition to the usual 
right-right and left-left bosonic correlators, there are also 
mixed ones, Eqs. (A9), (A19), and (A20). As a conse- 
quence there are four non vanishing 2-points correlators 
for ip a , instead of just two as for the full line. For in- 
stance, considering the critical case 77 = 0, when the U(l) 
transformation (23) is preserved, we have 

{lpl{h,Xi)lpi(t2,X 2 )) = (lpl(tl,Xi)lpl(t2,X2)) = 

[D(vt 12 - x 12 )f [D(vt 12 + x 12 )f 

x [V(vt 12 - x 12 )V{vt 12 + x l2 )] aT , (52) 

(^(il,£l)^2(*2,2;2)) = (lp2(tl,X 1 )lf>l(t2,X2)) = 

[v(vt 12 - Xl2 )r [v(vt 12 + Xl2 )r 

x [V(vt 12 - x 12 )f 2 [V{vt 12 + x 12 )f , (53) 

and 

(V>2 Zl)V>2(*2, X 2 )) = {^2{tl,X 1 )i)l{t2,X2)) (54) 

= (52) with a <-> t 

(^2(*l> X\)l\}\{t2, X 2 )) = {ll>l{tl,Xi)l/>2(t2>%2)) (55) 

= (53) with a <-» r , 

with x\2 = X\ + x 2 . The non-triviality of the correlators 
(53) and (55) reflects the breaking of the [/-symmetry on 
the half-line for 77 = 0. 

All the correlation functions just derived must be com- 
pared with the general scaling form coming from bound- 
ary conformal field theory 65 that in imaginary time Tj = 
iti predicts in general 



<**(*!)*(**)> = 
with the four point ratio 



1 



212^12 



(fan 



^11^22 
Z V1 Z 2\ 



F(0 , (56) 



(57) 



and Zi = Xi + irj, zi = Ei. F(t;) encodes all the boundary 
dependence and for small argument can be written as 65 
F(£ <C 1) oc £ db , where d b is called boundary exponent. 
The real time correlations we wrote are clearly not of this 
form, but this is just because we wrote them in the regime 
xi,x 2 S> 1 and x\ 2 , X12 arbitrary using the definitions 
(6) and (7) . If we want to get the correct scaling also for 
arbitrary xi. 2 we should modify the definitions as 

ipi{t,x) cx : e ^fR(vt-x) .. e i^r VL (vt+x) . ) (- 58 ) 
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at the price of introducing some more divergences that 
are easily renormalized. With this prescription, we ob- 
tain as a typical example 



{lpt(ti,Xi)lpi(t 2 ,X2)) = {lpl{tl,Xi)lp*(t2,X2)) 

[D(vt 12 - x X2 )f [D(vt 12 + x 12 )f 
T>{vt\2 — xi 2 )V(vt 12 + x\ 2 



V(2 Xl )V(2x 2 ) 



(60) 



that agrees with the general conformal field theory scal- 
ing with F(£) = ^°" r and so db — err. All the other corre- 
lation functions are easily modified accordingly. Because 
it will be easier to write, in the following we will ignore 
the double normal product and still use definitions (6) 
and (7). The expressions taking into account the cor- 
rect normalization at the boundary can be easily written 
down from the correlation we will derive. 

We finally point out that for Dirichlet boundary condi- 
tions, i.e. rj = oo, the diagonal correlations are the same 
but with df, = —<jt. Non diagonal correlations can be 
found in Ref. 34. 



C. Generic junction 

The case of a junction with an arbitrary number n > 1 
of wires can be actually reduced to the study of a suitable 
family of n half- lines. In fact, let U be the unitary matrix 
diagonalizing U which defines the boundary conditions 
(32). Since U is symmetric, we can choose U orthogonal, 
U l = U^ 1 , and real, U* = U. Let us parametrize the 
diagonal form 



V d =U Vit- 



as follows 



U d = diag (c- 2iQl ,e- 2iQ2 , 



-2ia T 



(61) 



) • (62) 



Using the definition (34) of S(k), one easily verifies that 
U diagonalizes S(k) for any k and that 



S d {k)=US(k)U- 1 



diag 



k — irji k — \r] 2 
k + ir/i ' k + \rj 2 



*0n 



«7r; 



where 



rn = Atan(ai) , 



7T 7T 
< Oi l < - . 

2 - 1 ~ 2 



(63) 



(64) 



Therefore S(fc) is a meromorphic function in the complex 
fc-plane, whose poles are different from and are all lo- 
cated on the imaginary axis. The condition (38) implies 
absence of bound states i.e. of poles in the upper complex 
fc-plane, namely < a-i < tt/2, hence rji > 0. 

Critical boundary conditions correspond to a matrix 
U such that the scattering matrix is insensitive to the 
momentum scale transformations A — > g\ (or k — > Q~ l k) 



with q > 0. To be scale invariant, the scattering matrix 
must have each rji vanishing or infinite, so that S is ac- 
tually momentum independent and with eigenvalues ±1. 
By means of Eqs. (35), (36), and the derivative 21 



_rfS(fc) 



1 



;(jfe)-s*(fe)]s(jfe) 



dk 2 

we see that criticality is equivalent to the condition 

§ = §* . 



(65) 



(66) 



In appendix B some examples of critical junctions with 
two, three and four wires are given. 

The matrix U allows us to define real scalar fields 
(p d = Uip which arc not localized on the single edges 
but have simple boundary conditions, formally the ones 
of disjoined half-lines 



0. 



l,...,n. (67) 



Comparing with the half- line Eqs. (49), it is straightfor- 
ward to derive the commutation relations for the right 
and left movers on the wires as done in Refs. 20-22 and 
reported in the appendix A. 



D. The TL model at the junction 

The TL model on the star graph T is defined by the 
sum of n Hamiltonians in Eq. (2) plus the boundary term 
that we implement through Eq. (32) at the bosonic level. 
The charges on each wire are defined via Eq. (3) and 
generate the U(l) ® U(l) phase transformations (A16) 
and (A17) leaving the Hamiltonian invariant. The cor- 
responding quantum equations of motion in the bulk are 
given by Eqs. (5) for each wire independently. 

In analogy with Eqs. (6) and (7), the solution of the 
equations of motion is given by the vertex operator 

^i (*,«,*) oc : e i ^WfiMvt-x)+r Vi , L (vt+^] :; ( 68 ) 

where the normalization constants are given in the ap- 
pendix A and depend on the anyon Klein factors. All 
bulk relations (the value of a, r and v, the form of the 
currents etc.) of TL model on the line are still valid for 
half infinite wires jointed in a single vertex. 

It is interesting to rewrite the boundary conditions (32) 
in terms of physical quantities of the model: in particular 
at the critical points (66) where <pn,{0 = <$V-l(£) (i- c - a 
generalized version of the unfolded picture of the half- 
line), the boundary conditions get a very simple form 



i±(t,0,») = T^S ij j±(t l ) i) 



(70) 



which simply fixes the splitting of the currents at the 
junction. Comparing this expression with the Kirchhoff 
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conditions (39) and (40), we see that at least one of two 
charges Q and Q, associated to p + and respectively, 
is dissipated by a non trivial flow at the vertex. Since p + 
generates the electric charge for the ip, (A16), we typi- 
cally require the Kirchhoff 's rule (39) to preserve electric 
charge, while Q conservation is lost. 

As for the half-line, the non trivial behavior of right- 
left correlators due to the presence of vertex, allows more 
non vanishing correlation functions with respect to the 
line case. Let us consider the two-points function for ip in 
the Fock representation, and let us focus for simplicity on 
the case of critical boundary conditions (66). Imposing 
the Kirchhoff 's rule on the charge Q generated by U(l), 
there are four non vanishing two-points correlators: 



tions 

(^t(t 1 ,x 1 ,i 1 )ip 1 (t 2 ,x 2 ,i 2 )) = ^^r- 

[V{vt 12 +x 12 )Y 2s ^ [V{vt 12 -x 12 )V{vt 12 +x 12 )Y T ^ , 

(75) 

(ipi(t 1 ,x 1 ,i 1 )ip 2 (t 2 ,x 2 ,i 2 )) = ZllZt2 

Ztt 

A [( ff 2 +r 2 )S lls2 +2 CT r5 nl2 ] [I3( ^ i2 _ Xi2)I3 ^ i2+Xi2) par5 llI2 

[D(vt 12 - x 12 )]-° 2 ^ [D(vt 12 + x 12 )]- t2s ^ , (76) 



(ipl(ti,xi,ii)ip 1 (t 2 ,x 2 ,i 2 )) = 



2ir 



[D(vt 12 +x 12 )f 5 ^ [V(vt 12 ~xi 2 )V{vt 12 +x 12 )Y T ^\ 



(71) 



{i>* {h,xi,ii)^ 2 (t 2 ,x 2 ,i 2 )) 



2ir 



K -[(° 2 +T 2 )^ 2+2aT s^ 2 ] [v{vti2 _ j 12 )]- 2 S ili2 

[V(vt u -x 12 )V(vt 12 +x 12 )Y rS ^[V(vt 12 +x 12 )Y 2s '^, 

(72) 

with all normalization factors defined in appendix A. All 
other non- vanishing correlation functions have the same 
form as the ones on the half- line Eqs. (52), (53), and 
(55) with only the proper wire index added. 
For the charge densities one finds 

(p + (t 1 ,x 1 ,i 1 )p + (t 2 ,x 2 ,i 2 )) = 
(2 ~^ )2 { [D 2 {vt 12 - x 12 )+V 2 (vt l2 + afia)] S lll2 

+ [V 2 (vt 12 -x 12 )+V 2 (vt 12 +x 12 )]§ lll2 } , (73) 
and for the currents 

(j+(h,x 1 ,i 1 )j+(t 2 ,x 2 ,i 2 )) = 

rz } s 9 ( \V 2 (vt 12 - x 12 ) + V 2 (vt 12 + x 12 )] S ill2 
(2ttC+) l 

- [V 2 (vt 12 -x l2 ) +V 2 (vt 12 +x 12 )] Si 1<a }. (74) 

The opposite signs in the Si 1 i 2 and Sj^ contributions in 
(74) ensure the Kirchhoff's rule for Q. Analogous expres- 
sions hold for p_ and j_ up to replace in Eqs. (73) and 
(74) (r + a) <-» (r - a) and § <-> -S. 

If instead we impose the conservation of the charge 
Q we have the non- vanishing two-point correlation func- 



(■4>i(t 1 ,x 1 ,i 1 )ipt(t 2 ,x 2 ,i 2 )) = 

(ipl(t 1 ,x 1 ,i 1 )il) 1 (t 2 ,x 2 ,i 2 )) , 
(4> 2 {t 1 ,xi,ii)ipi(t 2 ,x 2 ,i 2 )) = 

(tp 2 (ti,xi,h)ipl(t 2 ,x 2 ,i 2 )) = 
(ipl(ti,xi,h)tp 2 (t 2 ,x 2 ,i 2 )) = 

(ipi{t 1 ,x 1 ,i 1 )ip 2 (t 2 ,x 27 i 2 )) , (77) 

and 

(i(j 2 (t 1 ,x 1 ,i 1 )ip 2 (t 2 ,x 2 ,i 2 )) = 
{il) 2 {ti,xi,ii)il) 2 {t 2 ,x 2 ,i 2 )) = (75) with a <-> r . (78) 

The non conservation of the electrical charge is explicitly 
shown by the presence of non-neutral correlator (tpip). 
The correlations for conserved density p_ and current j_ 
are the same as Eqs. (73) and (74). 

IV. RG FLOW ON THE JUNCTION 

We completely characterized the fixed-point structure 
for a junction with an arbitrary number of wires n. Let us 
recall the main features explained in the previous section 
and in the appendix B. At the critical point, the scatter- 
ing matrix can only have eigenvalues ±1. For generic n, 
the fixed points are classified in terms of the integer num- 
ber p with < p < n, which is the number of eigenvalues 
equal to —1. At the fixed point, the boundary couplings 
rji (with 1 < i < n) are zero if the corresponding eigen- 
value is +1, infinity if the eigenvalue is —1. p = cor- 
responds to Neumann boundary conditions on all wires, 
while p — n to Dirichlct. Other values of p correspond 
to intermediate boundary conditions, that are n — p Neu- 
mann and p Dirichlet fields in the basis (pf diagonalizing 
the S-matrix. In Fig. 2 we report as a typical example 
the RG flow diagram for three wires in the r\i space. The 
final point of any axis is 77^ = 00. Let us discuss now 
the structure of the fixed points, postponing the study 
of the stability to the following. There are 2 3 = 8 fixed 
points families, one Neumann, 3 points with p = 1, three 
with p = 2 and one Dirichlet [in the general case, there 
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are 2™ families of which (™) for any p}. Every critical 
point belongs to a continuous family with p(n — p) real 
parameters that arc not shown in Fig. 2. Summarizing 
any critical point is identified by p, by the specific eigen- 
values that are —1 (i.e. by the axis in the figure) and by 
the p(n — p) real parameters. The parameters specifying 
the fixed point in the families are the angles cti reported 
for some examples in appendix B. For a given situation, 
the fixed point value of a, is given by their initial values. 
This means that ai are marginal couplings and their val- 
ues cannot be fixed only by requiring scale invariance. 

The role played by the conservation rules in this flow 
diagram is fundamental. To consider the most physi- 
cal case, let us discuss when the electrical charge is con- 
served, i.e. the Kirchhoff rule X^J+(A0,«) = is satis- 
fied. The first effect is to fix to zero one (arbitrary) 
constraining the system on the shadow area in Fig. 2 
so that Dirichlet boundary conditions are ruled out for 
the problem. Also the number of real parameters char- 
acterizing the p-fixed points is largely reduced. For three 
wires, the point with p = 1 becomes a one-parameter 
family, while the point with p — 2 becomes an isolated 
fixed point. Details for the general case are in the ap- 
pendix. Needless to say that imposing the conservation 
of the U(l) charge, results in fixing one of the r\i to oo and 
similarly reduced the number of real parameters available 
for each fixed point. 

We briefly discuss our terminology for the fixed points, 
in order to make the comparison with other papers as 
simple as possible. The fixed points with p = 2 ("2" in 
Fig. 2) is the mixed fixed point found by Nayak et. al 1 
and called D (or Dp) in Ref. 9 because of the n — 1 
Dirichlet boundary conditions (there n = 3) on the neu- 
tral modes (but this point is obviously different from our 
D). The family with p = 1 in Fig. 2 depends on a contin- 
uous real parameter a, as shown in Eq. (B4), and it has 
been first found in Ref. 21. Note that it is not symmet- 
rical under wire permutations. There are three special 
values of a: for a = — 1, 0, oo the S*-matrix breaks into a 
lxl and a 2 x 2 blocks. The lxl block is a wire decoupled 
from the other two that form a purely transmitting n = 2 
junction (the same can be verified for higher ra, changing 
the a's we can decouple any wire). For these special val- 
ues of a, the fixed points were also found by Chamon et 
al. 9 that called them asymmetrical Da- Other values of a 
interpolate continuously between these three. Finally it 
is worth commenting that the Dirichlet fixed point (D in 
Fig. 2) physically corresponds to n wires with an end in- 
serted into a large superconductor. In fact, the S-matrix 
S = — I gives conductance G = 21 corresponding to An- 
dreev reflection in all wires (i.e. sending a particle one 
gets an hole out). This is a different problem from a junc- 
tion of wires (even superconducting), because the large 
superconductor breaks the U(l) charge conservation 9,25 . 

Now we know the fixed-point structure, but what is 
the relative stability? Which fixed point describe the 
universal low energy behavior? There are several equiva- 
lent ways to tackle this question. The more natural one, 



as done elsewhere 1,9,23 , relies on calculating the scaling 
dimension of the perturbing operator at a given fixed 
point. Since our problem can be thought as n inde- 
pendent half-lines with n — p Neumann boundary con- 
ditions and p Dirichlet ones, the problem is just equiva- 
lent to understand the stability of Neumann or Dirichlet 
against a Robin term as in the Hamiltonian (46). This 
is a standard problem. In the bosonic theory, the flow 
can be followed exactly from Eq. (65) of the off-critical 
S'-matrix. The Neumann fixed point is always unstable, 
while Dirichlet is stable (or mixed if Kirchhoff is imposed 
on the electrical charge). However, as well known, con- 
sidering the fermionic theory changes this scenario be- 
cause of the Klein factors. In boundary conformal field 
theory, the stability conditions are just read from the 
boundary dimensions db appearing in the two-points cor- 
relation functions reported above. At the Neumann BC 
we have that the dimension is <jt = (£ 2 — C-)/4 that is 
greater than zero for g + > g_, i.e. for repulsive anyonic 
interaction, giving a stable Neumann. Oppositely at the 
Dirichlet BC the boundary dimension is —err that it is 
stable in the complementary attractive case. Since there 
are no other fixed points in the RG diagrams, this anal- 
ysis fixes all the RG flow. Note that for free anyons (and 
in particular fermions) r\ is marginal in this approach. 
In any given anyonic/fermionic model the actual stable 
fixed point will be determined by the higher order terms 
in r\ neglected in our approach. 

These results can be confirmed on the basis of the fol- 
lowing argument based on the so called ^-theorem 67,68 . 
For a one-dimensional critical system with a boundary, 
it is known that the boundary contribution to the en- 
tropy In g (g is the so called "universal non integer ground 
state degeneracy" 67 ) decreases along the renormalization 
group flow. We can easily calculate the value of the 
effective-potential V^ff = g+p+ + g-p- for the off critical 
model for any r\. Subtracting the divergent contribution 
of the bulk to make this expectation value finite, we get 
on each wire 

e(x,i) = (V eS (t,x,i)) = fi / — \k\e 2lkx S u (k) (79) 
where 

fully encodes the bulk interactions effect. In particular, 
when g + = g_ it vanishes and changes sign, giving the 
correct stability scenario. 

In fact, we can rewrite (79) in terms of the potential 
e Vj (x) for disjointed half-line with the boundary condi- 
tion (67) 

n 

£(*,*) = ]Ti^i%», (si) 

3 

with 

e v (x) = - A [1 - 4{x V ) - 8(x7 7 ) 2 e 2a; ' ) Ei(-2x? ? )] . (82) 
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FIG. 2: (Color online) RG flow diagram for a junction of 
three wires in the 771,2,3 space. The fixed points are: D is 
Dirichlet and corresponds to all 77; = 00; N is Neumann with 
r\i — 0; 1 are three fixed points families (depending on two 
parameters) with two 77 vanishing and one infinite; 2 are three 
fixed points two-parameter families with one 77 zero and two 
infinite. The cyan-shaded area is the allowed region when the 
Kirchhoff's rule for the electric charge is valid. It includes 
Neumann, two p = 1 families (with only one parameter left 
free) and one p = 2 fixed point (with no free parameter left). 
The arrow in the flow corresponds to the attractive case with 
g > 1 which gives Dirichlet as the most stable fixed point 
(without Kirchhoff) or the mixed p = 2 (with Kirchhoff). In 
the opposite repulsive case g < 1, all the arrows are reversed 
and the most stable fixed point is Neumann. 

The function 

n n 

s(x) = -Ax 2 i) = -4x 2 ^ e Vj (x) , (83) 

i=i 3=1 

collects the contribution of all the wires. It is a 
monotonous function with fixed points at 77 = 0, 00 in 
agreement with the g theorem. The stability of the fixed 
points and the direction of the flow are just given by the 
sign of f2 and agrees with the previous analysis. 



V. CONCLUSIONS 

In this paper we presented a systematic study of the 
critical properties of n anyonic Luttinger wires jointed in 
a single vertex. Imposing the boundary conditions (32) 
at the junction directly on bosonized fields allowed us to 
describe completely the RG flow diagram for any n. As a 
typical example the RG flow for n — 3 is depicted in Fig. 



2 where the main features of the various fixed points are 
discussed in the text. 

At this point it is worth comparing our findings with 
the literature. For two wires, our results are a simple 
anyonic generalization of the well-known ones by Kane 
and Fisher 29 for fermions that are reproduced for k = 1. 
For n — 3, as we said in the introduction the literature is 
enormous. The boundary conditions we used are equiv- 
alent to those of the "auxiliary model" of Nayak et al. 
for g 7^ 1 [in fact, expanding the exponential defining the 
auxiliary model 1 and keeping only up to the quadratic 
terms, neglecting irrelevant higher orders, we arrive to 
the Hamiltonian (46) where the symmetry of the bound- 
ary terms is just the Kirchhoff's rule]. We predict two 
possible stable fixed points: Neumann and mixed. Neu- 
mann is well known, it has zero conductance and in this 
setting it is stable for all repulsive interactions, i.e. g < 1. 
The mixed fixed point has been found for the first time 
by Nayak et al. 1 and it is specific of the junctions. It has 
enhanced conductance G/Gu nc = 4/3 and we found it is 
stable for all attractive interactions g > 1 as in Ref. 1. 
Everything agrees with the auxiliary model, but not with 
the "standard model" defined by the boundary condition 
(4), that is known to be different 1 . In fact in the standard 
model, the Neumann fixed point is stable only for g < 1/3 
while the mixed one only for g > 9. In the other regimes 
with 1/3 < g < 9, new fixed points appear that cannotbe 
present in our approach 1,9 . Our setting however presents 
a great advantage: it is simpler for generic n and more 
efficient in describing the off-critical properties of the sys- 
tem. In fact we provide for the first time the critical be- 
havior for all n. We found for g < 1 a Neumann stable 
fixed point (with zero conductance) and for g > 1 a mixed 
fixed point with conductance G/Gu nc — 2(n — l)/n. We 
also find other fixed points (described in the appendix B) 
that however have at least one direction of instability in 
the rji space and so they are multicritical points, in the 
sense that some other constraints must be imposed to 
reach them. Clearly we expect that the standard model 
for n > 3 will have some fixed points not found here, as 
for the case n — 3. A part from the per se interest of 
the model, the fixed points we found are relevant for the 
standard model as well. In fact, it is easy to generalize to 
any n the strong and weak boundary coupling (i.e. our 77) 
calculations of Refs. 1,9 to show that for small enough g 
the relevant fixed point is Neumann and for large enough 
is the mixed one. However, which fixed point governs the 
dynamics when none of these two is stable is not accessi- 
ble to our approach. For n = 4, two fixed points derived 
in the Appendix B have been recently found to describe 
the scattering matrix for a proposed experiment to de- 
tect the helical nature of the edge states in quantum Hall 
systems 26 . 

We mention that we also characterized the junction in 
the absence of the Kirchhoff's rule for the electric charge. 
It is of particular relevance considering the case when re- 
laxing the conservation of the electrical charge and im- 
posing the conservation of the dual one U(l). In this 
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case the more stable fixed point is always Dirichlet with 
uncommon points like the mixed one representing multi- 
critical points. 

There are two generalizations of the model considered 
here that should be easily accessible to a similar analysis. 
First of all one can consider fermions with spin (and even 
multispecies anyons) as done elsewhere with fermionic 
boundary conditions 23 . In this way one can understand 
which fixed points are present also with bosonic bound- 
ary conditions. The other generalization is relaxing the 
symmetry for time reversal to allow a non vanishing flux 
at the junction 9 . 

We close this paper on a more speculative level. In re- 
cent times there has been an increasing interest in quan- 
tifying the entanglement in extended quantum systems 
(see e.g. 69 as reviews). Among the various measures, 
the so-called entanglement entropy has by far been the 
most studied. By partitioning an extended quantum sys- 
tem into two blocks, the entanglement entropy is defined 
as the von Neumann entropy of the reduced density ma- 
trix pa of one of the two blocks. This procedure requires 
an arbitrary division of the system in two parts. In the 
junction problem studied here the system is automati- 
cally divided in parts and it would be very interesting 
to understand the amount of entanglement between the 
various wires. The analysis of some models on the line 
with one defect 70 (i.e. n = 2 in the language of this pa- 
per) showed that the entanglement entropy is not always 
only dependent on the central charge of the bulk the- 
ory (as maybe naively expected). The natural question 
is whether the conformal field theory formalism that has 
been successfully applied to the bulk and boundary case 71 
can be generalized to the junction. Furthermore, if we 
would be able to solve the non-equilibrium problem with 
changing the boundary condition (e.g. suddenly adding 
or removing the junctions, as done for n = 2 in Ref. 72), 
one can think of using the junction as an entanglement 
meter following the recent proposal based on quantum 
noise measurement 73 . 
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APPENDIX A: BOSONIZATION AND 
QUANTIZATION OF THE TL MODEL 

1. The line 

The basic tool for quantizing the system, described by 
the Eqs. (5), is the algebra A generated by the bosonic 



annihilation a(k) and creation a*(k) operators satisfying 

[a(k),a(p)] = [o*(fc),o*(p)] = 0, (Al) 
[a(k) , o*(p)] = 47r|fc- 1 | A< 5(fc - p) , (A2) 



where the normalization can be fixed such that 

feeT f 



\k- L \ A = ^ 



d_ 

dfc 



6{k) In- 



(A3) 



The derivative here is understood in the sense of dis- 
tributions, 7b is Euler's constant and A > is a free 
parameter with dimension of mass having a well-known 
infrared origin. It is useful to introduce 



u(A£) 



o 



AC 



-ifc£ 



-Iln(A|C|)-i e (0 

7T A 



= ln(iA£ + e) , e > . (A4) 

The left and right chiral fields are defined by 
C°° c\k 

MO = J o —[a*(k)e ik t+a(k)e- ik t] 
r°° Ah. 

MO = J —[a*(-k)e^+a(-k)e 
and obey the commutation relations 
[Mb) , Mb)} = [Mb) , Mb)} 



-ifcfi 



(A5) 
(A6) 



-i£(& 2 ), (A7) 



[Mb),Mb)} = 0, (A8) 
and have the correlations 

(Mb)Mb)) - (Mb)Mb)) = u(A& 2 ) , (A9) 

with £12 = £i - £ 2 and obviously (<^fl(£i Vl(6)) = 0. 
Defining the chiral charges by 

Qz = \l d£(0¥>z)(O. Z = R,L, (A10) 
one gets 

[Qr,M0} = [Qi.VL(0] = -V2, 
[Qr,M0] = [Ql,MO) = IQr,Ql] = o.(aii) 

It is worth mentioning that all previous the commutation 
relations are invariant under the duality transformation 

MO -> MO . MO -> -MO . (A12) 

which define the T-duality in string theory. 

At this point we are ready to introduce a family of 
vertex operators parametrized by two real variables a 
and r defined by 

A(t,x) — ze l ^( T Q R ~ a Q L ' > • c i Vn[c<PR(vt~x)+T<p L (vt+x)] . 

(A13) 
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with 



z = (2^)- 1 /2 A (- 2 +- 2 )/ 2 : 



(A14) 



where : • • • : denotes the normal product in A and v is 
some velocity to be determined by consistency. From 
Eqs. (6) and (7) the fields ipi and ip 2 are vertex oper- 
ators with interchanged a and r, with a normalization 
constant given by Eq. (A14). The factor ^(rQ R -<^Q L ) 
is included in the definition (A14) to ensure canonical 
anionic commutation relation between ipi i2 without in- 
troducing Klein factors that will be important only for 
the fields on different wires. 

The following identity is useful in determining the ex- 
change properties of the vertex operators and so all cor- 
relation functions 



A*(t,x 1 )A(t,x 2 ) = \x 12 \ 



>e 2' 



: e 



iy/TT[a^p R (v t — X 2 ) — Glfi R (vt — Xt) + T<p L (vt + X 2 ) — TifiLivt + Xx)} . 



(A15) 



where x\ 2 = x\ — x 2 . 

The normalization of the charge densities p± is fixed 
by requiring that they generate the transformations (24) 
and (23) in infinitesimal form, namely 

[p+(t,xi),ip a (t,x 2 )]= -${x 12 )ip a (t,x 2 ) , (A16) 

\p-{t,x 1 ),%l) a {t,x 2 )\= -{-l) a 6{x 12 )^ a {t,x 2 ). (A17) 

2. The half-line 

In the main text, we stressed that on the half line right 
and left modes couple and have non trivial commutation 
relations given by Eq. (49). This gives rise to few changes 
to the relations valid on the full line. The vertex operator 
is always defined by Eq. (A13), but the normalization 
constant is affected by the boundary 34 : 



(27r)- 1 /2 A (-+-) 2 /2 ) ^ = 0- 



r (27T)- 1 
1 (27T)- 1 



/2 A (a-r) 2 /2 5 < 7? < OO . 



(A18) 



The right- left coupling also affects the correlation func- 
tions of the field (p. In fact, while the right-right and 
left-left correlators are still given by Eq. (A9), the mixed 
ones are 



(v?-r(£iVl(&)) = 



u(A£i 2 ) r) = 0, 

-u(A£i 2 ) r] = oo, 

-u(A£ 12 ) - V-(A£ 12 ) < rj < oo, 

(A19) 



u(A£i 2 ) r) = 0, 

(<Pl(£i)<Pr{&)) = { -u(Hi2) 77 = 00, 

-u(A^i 2 ) - v+(A£i 2 ) < r) < oo, 

(A20) 

where the "boundary propagator" is 



IT 



(A21) 



and Ei(x) = dze~ z / z is the exponential integral func- 
tion, that at small x has the right logarithm expansion. 
Note that in the above formulas for mixed correlators 
£i = vt\ — x\ and £ 2 = vt 2 + x 2 or viceversa, thus 
£12 = vt\ 2 =F x\ 2 , with the sign depending on the cor- 
relator if it is right-left or left-right respectively. 



3. The junction 

For the theory on the star graph, all the relevant com- 
mutation relations and correlators of the fields follow 
from those on the half line after performing the linear 
transformation U in Eq. (61). In fact all the fields <p d 
are just delocalized fields satisfying the proper bound- 
ary conditions reported above with different rji for each 
mode. Thus, comparing with the half line equations (49), 
it is straightforward to derive the commutation relations 
for the right and left movers on the wires 

bi!,fi(£i), ^2,^(62)] = [<P»i,l(£i),<P*2 ,£(&)] 

= -^12)*^, (A22) 

[PnAZ^ViiAb)] =^ 1 ^ 2 [<.i?(a),<,L(6)] , 

(A23) 

where v4 )£ (0 = U<Pr,l{Q and 

[<*(£i),<l(6)] = 

-ie(£i2)£iii 2 , Vh = 0; 

H^)Si li2 , % = 00 ; 

[ie(&2) - 4i0(a 2 )e-' ) - 1 « 12 ] 5 lll2 , < Vil < 00 . 

(A24) 

The mixed commutator (A23) simplifies greatly for crit- 
ical boundary conditions 

Vi a ,L(6,*2)] = -ie(£i 2 )§ui 2 • (A25) 

Note that at spacelike distances where vt\ 2 — x\ 2 < 0, the 
commutators (A22) and (A23) behave as if the scattering 
matrix were replaced by the critical one obtained in the 
infrared limit A — > 00 or equivalcntly k — > 

K,ii(£i),W 2 ,L(£2)]|, 2t2 x2 =-iSi li2 (0). (A26) 

This simply means that tpn.L has the same properties of 
locality than its infrared limit. 

The last complication on the star graph arises in the 
definition of the anyonic fields ipi, 2 - To have the correct 
commutation relation they must be defined according to 

V>l(t, X, i) = Z, Tli jV^irQi.R-crQi.L) 

X • e iVK[a<Pi,R(vt-x)+Tipi iL (vt+x)] . 

ih(t, x, i) = z t m e i v^( CT< 5>,H- 1 -<3^) 

x : e i V^l T v>i,R(vt-x)+crip iiL (vt+x)] . ^ (A27) 
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where Zi are fixed to 

z t = (2 7 r)- 1 /2A[( cr2 + T2 )+ 2 ^ s "(°)]/ 2 . 



(A28) 



and rji are the anyonic Klein factors needed to ensure the 
correct commutation of anyon fields on different edges 



ip(t,Xi,i)ip(t,Xj,j) = e 



(A29) 



parameters a\_ 2 : 



S 2 («i, a 2 ) = 



1 



1 + a{ + aj. 



a\ — a\ — 1 2a\a 2 2ct\ 

2aia 2 -al + al-l 2a 2 |, (B3) 
2ai 2a 2 1 - a\ - a\ 



where = e(i — j). ft is straightforward to build them and 
for example in terms of the auxiliary Majorana algebra 
[cj, Cj] = ineij and c* = q resulting in r\i =: c 7 "^ :. These 
factors are of fundamental importance when considering 
as junction condition Eq. (4), because it is written in 
terms of anyonic degrees of freedom. Oppositely, because 
the junction condition we use is written in terms of cur- 
rents that only get (re)normalized by the statistics, they 
are inessential. For this reason we do not discuss them 
further, remanding the interested readers to the complete 
treatment presented in Ref. 74 and in the appendix E of 
9. 



(B4) 



For generic values of the parameters these ^-matrices 
violate both U(l) and t/(l). Preserving t/(l), one must 
impose (39) on (B3). This implies ot\ = a 2 = 1, leading 
to the isolated critical point 



-12 2 
2 = ^1 2 -1 2 
6 1 2 2-1 



(B5) 



APPENDIX B: CRITICAL POINTS 

By scale invariance any critical point is associated with 
a k- independent S- matrix satisfying unitarity (35), Hcr- 
mitian analyticity (36) and time-reversal invariance (37), 
i.e. 



5* =S. 



3* - ! 



(Bl) 



The classification of these S-matrices is now a simple 
matter. Indeed, one can easily deduce from (Bl) that the 
eigenvalues of § are ±1. Let us denote by p the number 
of eigenvalues —1. The values p = and p — n corre- 
spond to the familiar Neumann (Sjv = I) and Dirichlet 
(Sd = —I) boundary conditions respectively. A richer 
structure appears for < p < n. In that case the S- 
matrices satisfying (Bl) depend onp(n- p) > 1 parame- 
ters, giving raise to whole families of critical points 20,21 . 
Let us describe them explicitly for n — 2, 3, 4. 

The only possibility for n = 2 is p = 1, leading to the 
one-parameter family 39,40 



1 



1 



-2a 1 - a 2 



(B2) 



with a a real number. For a = — 1 one has full transmis- 
sion and no reflection, which corresponds to the theory on 
the whole line. This is an example of exceptional bound- 
ary conditions already mentioned 66 . It is only the only 
S matrix in the family satisfying Kirchhoff 's rule for the 
electric charge. Oppositely, a — 1 is the only matrix sat- 
isfying Kirchhoff 's rule for the U(l) charge, as predicted 
by duality. 

In the case n = 3 one has two possibilities: p = 2 
and p = 1. In both cases one has a family with two real 



which is invariant under edge permutations. From (B4) 
one obtains instead a 2 = — (1 + a\). Therefore, setting 
a = a 2 , one has in this case the one-parameter family of 
critical points 



Si = 



a + a 



—a a(a + 1) 
, , a(a + 1) a + 1 



a + 1 



1 + a 

—a 

a a(a + 1) 



(B6) 

which is not invariant under edge permutations for 
generic a. Summarizing, the critical points which respect 
17(1) are S = I 3 , (B5), and (B6). The matrix (B5) has 
been discovered by means of RG techniques by Nayak et 
al. 1 . The family (B6) appeared for the first time in 21. 

If one wants to preserve U(l), one must require (40). 
One is left therefore with § 3 = — 1 3 , 



1 



a + cr 



—a a(a + 1) 1 + a 
a(a + 1) a + 1 —a 
a + 1 —a a(a + 1) 



and 



Si=-i 




(B7) 



(B8) 



as predicted by duality. 

For n = 4 the general matrices satisfying all the con- 
straints (Bl) are too large to be reported here. Thus 
we only give the critical points for n = 4 satisfying the 
Kirchhoff's rule Eq. (39) for the electrical current [the 
analogous ones with the Kirchhoff's rule Eq. (39) are just 
— S because of duality]. Besides S = I4 corresponding 
to p = 0, one has: 
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(i) for p = 1 the §-matrix depends on two real param- 
eters a\ 2 an d results to be 



333 = 



1 

AT 
1 

AT 
1 



{a\ + a\ + a 2 + a x a 2 + a\) 
(1 + ai + a\ + a 2 + aia 2 ) , 



(1 + ai + a 2 + a x a 2 + a\) , 



1 

a7 



*12 



514 



Al 



a 2 . 



(ai + a 2 + a\a 2 ) , 



§13 = --7— Oil 



— (1 + ai + a 2 ) , 
Ai 



1 1 

§23 = 7 OL\OL 2 , §24 = -t— a 2 (l + «1 + a 2 ) , 
Ai Ai 

§34 = -t— ai(l + ai +a 2 ) , 
Ai 

with Ai = 1 + ai + af + a 2 + a x a 2 + a 2 . The remaining 

entries are recovered by symmetry. Note that this matrix 
is not invariant under edge permutations. 

(ii) for p = 2 the S'-matrix still depends on two real 
parameters: 



§11 = — [3a 2 1 + 2a 1 (l-a 2 )-(l + a 2 ) 2 }, 

§22 = -^[-l-a? + 2a 2 + 3a|-2a!i(l + a!2)] 
A 2 

§33 - ~-[i-a 2 1 +2a 2 -al + 2a 1 {l + a 2 )], 

§44 = --^[a? + 2ai(l-a 2 ) + (l + a 2 ) 2 ], 
A 2 



51 2 = — — (1 + ai + a 2 + 2a\a 2 ) , 

A2 
2 

51 3 = — [a 2 (l + a 2 ) - ai(2 + a 2 )] , 

A2 
2 

514 = — — (1 + ai — aia2 + ot 2 ) , 

A 2 

323 = -r— (ai + ai - 2a 2 - aia 2 ) , 

A 2 

2 

3 2 4 = -7— (1 + ct\ + ol 2 - a x a 2 ) , 

A2 
2 

534 = -7— (ai + ol\ + a 2 + a 2 ) , 
A 2 



where A 2 = 3 + 3af + 2ai(l - a 2 ) + 2a 2 + Also 
this matrix is not invariant under edge permutations. 



(iii) for p = 3 the we have only an isolated S'-matrix: 



(B9) 



/ -2 2 2 2 
12-22 2 
4 2 2 -2 2 
V 2 2 2 -2 
which is invariant under edge permutation. This is the 
analogous for four wires of the Nayak et al. result 1 . 



Recently, the p = 1 matrix with ot\ = 1 and a 2 = — 1 
and the p = 3 matrix has been found to describe the 
scattering matrix for a proposed experiment to detect 
the helical nature of the edge states in quantum Hall 
systems. 



We conclude this appendix with the matrix with p = 
n — 1 for general n satisfying the electric Kirchhoff rule: 



1 

n 



I (2 - n) 2 2 ■■■ 2 
2 (2-n) 2 ■■■ 2 

v 2 2 2 ••• (2-n) J 



(BIO) 



since it is the most stable in the RG phase diagram as 
shown in the text. This is the only matrix which is invari- 
ant under wire permutations (i.e. that has all diagonal 
elements equal and non-diagonal as well), satisfying the 
Kirchhoff 's rule and with all non-vanishing entries. 
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